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Abstract 

Fast and efficient numerical-analytical approach is proposed for modeling complex collec- 
tive behaviour in accelerator/plasma physics models based on BBGKY hierarchy of kinetic 
equations. Our calculations are based on variational and multiresolution approaches in the 
bases of polynomial tensor algebras of generalized coherent states/wavelets. We construct 
the representation for hierarchy of reduced distribution functions via the multiscale decom- 
position in high-localized eigenmodes. Numerical modeling shows the creation of different 
internal coherent structures from localized modes, which are related to stable/unstable type 
of behaviour and corresponding pattern (waveletons) formation. 
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Abstract 

Fast and efficient numerical-analytical approach is pro- 
posed for modeling complex collective behaviour in ac- 
celerator/plasma physics models based on BBGKY hier- 
archy of kinetic equations. Our calculations are based 
on variational and multiresolution approaches in the bases 
of polynomial tensor algebras of generalized coherent 
states/wavelets. We construct the representation for hier- 
archy of reduced distribution functions via the multiscale 
decomposition in high-localized eigenmodes. Numerical 
modeling shows the creation of different internal coherent 
structures from localized modes, which are related to sta- 
ble/unstable type of behaviour and corresponding pattern 
(waveletons) formation. 

1 INTRODUCTION 

The kinetic theory describes a lot of phenomena in 
beam/plasma physics which cannot be understood on the 
thermodynamic or/and fluid models level. We mean first of 
all (local) fluctuations from equilibrium state and collec- 
tive/relaxation phenomena. It is well-known that only ki- 
netic approach can describe Landau damping, intra-beam 
scattering, while Schottky noise and associated cooling 
technique depend on the understanding of spectrum of lo- 
cal fluctuations of the beam charge density [1]. In this 
paper we consider the applications of a new numerical- 
analytical technique based on wavelet analysis approach 
for calculations related to description of complex collec- 
tive behaviour in the framework of general BBGKY hierar- 
chy. The rational type of nonlinearities allows us to use our 
results from [2]-[15], which are based on the application 
of wavelet analysis technique and variational formulation 
of initial nonlinear problems. Wavelet analysis is a set of 
mathematical methods which give us a possibility to work 
with well-localized bases in functional spaces and provide 
maximum sparse forms for the general type of operators 
(differential, integral, pseudodifferential) in such bases. It 
provides the best possible rates of convergence and mini- 
mal complexity of algorithms inside and as a result saves 
CPU time and HDD space. In part 2 set-up for kinetic 
BBGKY hierarchy is described. In part 3 we present ex- 
plicit analytical construction for solutions of hierarchy of 
equations from part 2, which is based on tensor algebra ex- 
tensions of multiresolution representation and variational 
formulation. We give explicit representation for hierarchy 
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of n-particle reduced distribution functions in the base of 
high-localized generalized coherent (regarding underlying 
affine group) states given by polynomial tensor algebra of 
wavelets, which takes into account contributions from all 
underlying hidden multiscales from the coarsest scale of 
resolution to the finest one to provide full information about 
stochastic dynamical process. So, our approach resembles 
Bogolubov and related approaches but we don't use any 
perturbation technique (like virial expansion) or lineariza- 
tion procedures. Numerical modeling shows the creation 
of different internal (coherent) structures from localized 
modes, which are related to stable (equilibrium) or unstable 
type of behaviour and corresponding pattern (waveletons) 
formation. 

2 BBGKY HIERARCHY 

Let M be the phase space of ensemble of N particles 
(dmxM = 6 AO with coordinates X{ = (qi,Pi), i = 
1,...,N, qi = (ql,qi, qf) G i? 3 , Pi = (j>l,Pi,Pi) & 
R 3 , q — (qi, . . . , g^v) G R 3N . Individual and collective 
measures are: 



N 



fii = dxi = dq t pi, M = JJ M> 



(1) 



i=l 



Distribution function Dn(x\, . . . ,XN',t) satisfies Liou- 
ville equation of motion for ensemble with Hamiltonian 

Hn ■ 



dD 



N 



Of 



— {Hn, Dn} 



(2) 



(3) 



and normalization constraint 

D N (xi, . . .,x N ;t)dfi = 1 
where Poisson brackets are: 

{^,^} = E(f^f^-f^f^) (4) 
^ V dqi dpi dp t dqi I 

Our constructions can be applied to the following general 
Hamiltonians: 



N n 
Pi 



H »=Y,(%t + u ^) + E u a(M) ( 5 ) 

<=1 l<i<j<N 



where potentials Ui(q) = Ui(q 1} . . . , q N ) and Uij{qi,qj) 
are not more than rational functions on coordinates. Let i, 



2 



and Lij be the Liouvillean operators (vector fields) 



_ dUa d dUij d 

J-iii — 



J%3 dqi dpi dqj dpj 



(7) 



For s=N we have the following representation for Liouvil- 
lean vector field 

L N = {H N , •} (8) 
and the corresponding ensemble equation of motion: 

dD N 



dt 



+ L N D N = 



Ln is self-adjoint operator regarding standard pairing on 
the set of phase space functions. Let 

F N (xx, . . .,x N ;t) = y^D N {xi, ...,x N ;t) (10) 

Sn 

be the N-particle distribution function (Sjv is permutation 
group of N elements). Then we have the hierarchy of re- 
duced distribution functions (V s is the corresponding nor- 
malized volume factor) 

F s (x 1 ,...,x s ;t)= (11) 
V s / D N (xi, ...,x N ;t) J] m 

J s+l<i<N 

After standard manipulations we arrived to BBGKY hier- 
archy [1]: 



dF s 
dt 



+ L S F S = - [ d/is+i "V] Li, s+1 F S+1 (12) 



It should be noted that we may apply our approach even to 
more general formulation than (12). Some particular case 
is considered in [16]. 

3 MULTISCALE ANALYSIS 

The infinite hierarchy of distribution functions satisfying 
system (12) in the thermodynamical limit is: 

F = {F ,F 1 (x 1 ;t),F 2 {x 1 ,x 2 ;t), . . . , (13) 
F N (x!, . . . ,x N ;t), . . .} 

where F p (xi, . . . ,x p ;t) G HP, H° = R, HP = 
L 2 (R 6p ) (or any different proper functional space), F G 
H°° = H° H 1 • • • ® H v © . . . with the natural Fock- 
space like norm (of course, we keep in mind the positivity 
of the full measure): 



(F,F)=F* + Y / f F?{x 1 ,...,x i ;t)f[to 

i J e=i 



(14) 



First of all we consider F = F(t) as function on time 
variable only, F G L 2 (R), via multiresolution decompo- 
sition which naturally and efficiently introduces the infinite 
sequence of underlying hidden scales into the game [17]. 
Because affine group of translations and dilations is inside 
the approach, this method resembles the action of a micro- 
scope. We have contribution to final result from each scale 
of resolution from the whole infinite scale of spaces. Let 
the closed subspace Vj (j G Z) correspond to level j of res- 
olution, or to scale j. We consider a multiresolution anal- 
ysis of L 2 (R) (of course, we may consider any different 
functional space) which is a sequence of increasing closed 
subspaces Vf ...V- 2 C V-i C V C Vi C V 2 C ... satis- 
fying the following properties: let Wj be the orthonormal 
complement of Vj with respect to Vj+i: Vj+i = Vj Wj 



then we have the following decomposition: 



{*"(*)}= Wj 



(15) 



— oo< 7 <oo 

or in case when Vq is the coarsest scale of resolution: 



{F(t)} = V ®W j , 

3=0 



(16) 



Subgroup of translations generates basis for fixed scale 
number: span feeZ {2^ 2 *(2^ - k)} = W r The whole ba- 
sis is generated by action of full affine group: 

span fceZjeZ {2^ 2 *(2^ - k)} = (17) 
span feijez {* j>fc } - {F(t)} 

Let sequence {Vj}, Vj C L 2 (R) correspond to multireso- 
lution analysis on time axis, {Vp } correspond to multires- 
olution analysis for coordinate x%, then 



V 



n+1 \rX 



Vf 1 <g> . . . <g) Vf n ® V* 



(18) 



corresponds to multiresolution analysis for n-particle dis- 
tribution fuction F n (x\, . . . , x n ; t). E.g., for n — 2: 

V 2 = {f:,f(x 1 ,x 2 )= (19) 
^2 a ki,k 2 4> 2 (xi - h,x 2 - k 2 ), a klM G £ 2 {Z 2 )}, 

ki ,k 2 

where 4> 2 {xi, x 2 ) = (j) 1 (xx)4> 2 (x 2 ) = (j) 1 <8> 4> 2 {xi, x 2 ), 
and (f> l (xi) = 4>(xi) form a multiresolution basis corre- 
sponding to {Vp}. If {^(x! - £)}, £ G Z form an or- 
thonormal set, then </> 2 (xi — k\, x 2 — k 2 ) form an orthonor- 
mal basis for Vq. Action of affine group provides us by 
multiresolution representation of L 2 (R 2 ). After introduc- 
ing detail spaces W 2 , we have, e.g. V 2 — V^ © Wfi . Then 
3-component basis for W§ is generated by translations of 
three functions 

* 2 = 4>\ Xl ) ® * 2 (x 2 ), * 2 = ^\ Xl ) ® 2 (x 2 ), 
* 2 = * 1 (x 1 )®* 2 (a; 2 ) (20) 



3 



Also, we may use the rectangle lattice of scales and one- 
dimentional wavelet decomposition : 

f(xi,x 2 )= < ®*j>* > Vj,t®yj,k(xi,X2) 

i,£;j,k 

where bases functions 9 a ® depend on two scales 
2~ l and 2~- J . After constructing multidimension bases we 
apply one of variational procedures from [2]-[16]. As a 
result the solution of equations (12) has the following mul- 
tiscale/multiresolution decomposition via nonlinear high- 
localized eigenmodes 

F(t,xi,x 2 ,...) = ^2 a v Ul ® V j (t,xi,x 2 , . . .) 



Vi(t)=Vfr 8l ° w (t) + J2Vi i (u>it), c^-2 1 (21) 

1>N 

U\x s ) = U% slow (x s ) + £ UHk^Xs), k s m ~ 2" 



m>M 



which corresponds to the full multiresolution expansion 
in all underlying time/space scales. Formulae (21) give 




Figure 1 : 6-eigenmodes representation. 

us expansion into the slow part ^^°m m ^ f ast oscillating 
parts for arbitrary N, M. So, we may move from coarse 
scales of resolution to the finest one for obtaining more de- 
tailed information about our dynamical process. The first 
terms in the RHS of formulae (21) correspond on the global 
level of function space decomposition to resolution space 
and the second ones to detail space. In this way we give 
contribution to our full solution from each scale of res- 
olution or each time/space scale or from each nonlinear 
eigenmode. It should be noted that such representations 
give the best possible localization properties in the corre- 
sponding (phase)space/time coordinates. In contrast with 
different approaches formulae (21) do not use perturbation 
technique or linearization procedures. Numerical calcula- 
tions are based on compactly supported wavelets and re- 
lated wavelet families and on evaluation of the accuracy 
regarding norm (14): 



\F 



1 -F N \ 



< e 



(22) 



Fig. 1 demonstrates 6-scale/eigenmodes (waveletons) con- 
struction for solution of equations like (12). So, by us- 



ing wavelet bases with their good (phase) space/time local- 
ization properties we can construct high-localized wavele- 
ton structures in spatially-extended stochastic systems with 
collective behaviour. 
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